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Abstract 


For the ASO-S/HXI payload, the accuracy of the flare reconstruction is reliant on important factors such as the 
alignment of the dual grating and the precise measurement of observation orientation. To guarantee optimal 
functionality of the instrument throughout its life cycle, the Solar Aspect System (SAS) is imperative to ensure that 
measurements are accurate and reliable. This is achieved by capturing the target motion and utilizing a physical 
model-based inversion algorithm. However, the SAS optical system’s inversion model is a typical ill-posed inverse 
problem due to its optical parameters, which results in small target sampling errors triggering unacceptable shifts in 
the solution. To enhance inversion accuracy and make it more robust against observation errors, we suggest 
dividing the inversion operation into two stages based on the SAS spot motion model. First, the as-rigid-as- 
possible (ARAP) transformation algorithm calculates the relative rotations and an intermediate variable between 
the substrates. Second, we solve an inversion linear equation for the relative translation of the substrates, the offset 
of the optical axes, and the observation orientation. To address the ill-posed challenge, the Tikhonov method 
grounded on the discrepancy criterion and the maximum a posteriori (MAP) method founded on the Bayesian 
framework are utilized. The simulation results exhibit that the ARAP method achieves a solution with a rotational 
error of roughly +3”5 (1 /2-quantile); both regularization techniques are successful in enhancing the stability of the 
solution, the variance of error in the MAP method is even smaller—it achieves a translational error of 
approximately +18 um (1/2-quantile) in comparison to the Tikhonov method’s error of around +24 um (1/2- 
quantile). Furthermore, the SAS practical application data indicates the method’s usability in this study. Lastly, this 
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paper discusses the intrinsic interconnections between the regularization methods. 
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1. Introduction 


With the ongoing advancements in technology, there has 
been a renewed enthusiasm for solar investigation. In recent 
years, China’s heliophysical research has advanced signifi- 
cantly, yielding numerous representative outcomes and execut- 
ing some space-based solar expedition missions with specific 
objectives (Gan et al. 2019a). The ASO-S mission was born 
under this opportunity, it is China’s inaugural synthesized solar 
observatory satellite focused on detecting magnetic fields, 
flares, and coronal mass ejections of the Sun (Gan et al. 2019b). 
HXI is one of the three primary loads designed specifically for 
detecting solar flares at hard X-ray spectra. Its imaging utilizes 
a space-modulated dual grating collimator instrument, as the 
principle of this imaging technique requires the exact relative 
positioning of the grating slits (Su et al. 2019) and the precise 
solar orientation simultaneously (Zhang et al. 2019), which can 
otherwise negatively impact the accuracy of the flare intensity 
distribution inversion. To ensure this, an online monitoring 


system, integrated with the instrument, is necessary to measure 
the collimator’s status and observation direction. The Solar 
flare detection payloads that have been launched are all 
contained comparable systems. For instance, the Solar Aspect 
System (SAS) in RHESSI achieves a solar pointing accuracy of 
0”4 while the Twist Monitoring System yields a torsion 
measurement accuracy of a few arcsec (Zehnder et al. 2003). 
The STIX payload, launched in 2020 onboard the Solar 
Orbiter, also featured an SAS with a solar pointing accuracy 
of +4” (Krucker et al. 2020; Warmuth et al. 2020). In HXI, for 
the sake of measurement convenience, we choose the mature 
visual measurement scheme, which has been widely applied to 
different fields and scenarios, such as object detection (Zou 
et al. 2023), parameters measurement (Hashmi et al. 2022), 
structure and deformation monitoring (Dong & Catbas 2021; 
Zhuang et al. 2022), 3D reconstruction (Ham et al. 2019), 
robots (Abdulazeez & Faizi 2021), and Autonomous Vehicles 
(Joel et al. 2020). The SAS system uses the images to deduce 
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the deformation and observation orientation of the collimator 
by tracking the movement of a light spot. However, the 
significant differences in the parameters of the two lenses of 
SAS lead to a minute eigenvalue of the coefficient matrix in the 
inversion function. This is a common ill-posed inverse problem 
(Kabanikhin 2008; Engl & Groetsch 2014), where slight target 
sampling errors may result in oscillations and non-convergence 
of inversion results (Kabanikhin 2008; Engl & Groetsch 2014), 
making it impossible to meet the accuracy requirements with 
conventional solution methods like Gaussian, Newtonian, or 
Least Squares Optimization. 

To address the aforementioned issue, two primary studies 
were conducted: first, utilizing the rigid body transformation 
model to mitigate solution instability caused by outliers; and 
second, implementing the regularization method to tackle the 
pathological inverse problem. To accomplish this objective, we 
commence with the physical model of SAS measurements and 
divide the solution process into two stages. Initially, we resolve 
the intermediate variables from a fraction of the light spots in 
the system. Subsequently, we solve the remaining variables 
together with other characteristic light spots. One advantage of 
using this approach is the utilization of the as-rigid-as-possible 
(ARAP) algorithm in the initial stage. This algorithm is a 
highly efficient 3D morphology editing tool that can ensure 
maximum maintenance of the local mesh shape characteristics 
and produce a solution closely resembling a rigid-body 
transformation (Sorkine & Alexa 2007). The equations 
resolved during the second step following decomposition 
undergo formal simplification. However, the unresolved ill- 
posed problem presents an obstacle, which curbs the precision 
of the conclusive resolution. To address this issue, regulariza- 
tion theory is employed, a prevalent approach for dealing with 
ill-posed problems (Poggio et al. 1987; Hansen 2010). This 
method involves appending a penalty term to minimize the 
residuals norm, limiting the solution’s scope so that it is no 
longer an “exact solution” but a harmonization of the residuals 
and the solution’s range. Although this solution is not the 
precise value when the equations are solved rigorously 
(assuming that the equations are linear and of full rank), this 
resolution is still more acceptable than the non-convergence of 
the solution resulting from the ill-posedness of the equations. 

To verify the method’s effectiveness, we make simulations 
for the algorithm mentioned above. The results indicate that it 
solved the ill-posed problem that we encountered to a large 
extent. Besides, the method is of significant practical value as it 
is implemented in the HXI instrument’s actual production 
process. It helps in collimator assembly and grating adjustment 
and also achieves inflight monitoring and _ solar-pointing 
functions. Our research has assisted SAS in overcoming the 
limitations of its principle and improving its functionality, 
ensuring the success of the HXI observation mission. 
Furthermore, our findings highlight the optimization benefits 
of the regularization theory for ill-posed problems. This 
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extends the potential application of visual measurement 
methods and provides valuable insights for the development 
of the next generation of solar observation equipment. 


2. Measurement Principle 


The SAS functions as a visual measurement system that 
solves the state of the HXI collimator through inverse analysis 
of light spot images. The system is comprised of two sub- 
optical systems, namely the deformation monitoring subsystem 
(DM) and the solar aspect subsystem (SA). The DM is a short- 
focus optical system that possesses a focal length of 32.8 mm, 
with a diagonal field of view (FOV) of 27°8. It captures images 
of both the three frosted glasses affixed to the front grating 
substrate and the Sun, producing four spots on the detector. The 
SA is a long-focus optical apparatus with a focal length of 
about 1214 mm and a diagonal FOV of 45’, exclusively capable 
of imaging the Sun. Contrarily, the DM varies from the SA in 
that the DM lens is mounted on the SAS electronic control box, 
integrated with the detector, while the SA lens is suspended on 
the front substrate and positioned apart from the detector. The 
mounting configuration facilitates monitoring of the SA’s 
optical lens by the DM as a feature point, thus constructing an 
additional equation to create a more complete physical model 
that represents a more realistic state of the collimator. Figure 1 
displays the mounting state of SAS in HXI. Figure 2 explains 
how SAS works. 

We identify four variables that need to be concerned when 
monitoring the collimator state, which: the relative rotation a 
of the front and rear grating substrates along the collimator 
main axis; the relative translation D perpendicular to the main 
axis between the front and rear grating substrates; the tilt offset 
b of the DM optical axis; and the tilt offset y of the SA optical 
axis (also calibrated as the main axis of the collimator). We 
establish the motion pattern of the light spots to determine the 
four variables, of which a and d have a direct influence on the 
imaging effect on flares of the collimator and therefore are the 
main variables to be monitored by the SAS system. While 3 
and y are also the state variables of the HXI, but require less 
necessity for monitoring compared to a and D. It is mainly due 
to the fact that 8 represents the tilt of the SAS’s optical axis, 
which does not influence the imaging effect, and the y, while it 
affects the flare reconstruction by representing the Sun’s 
pointing direction, its value can be attained through alternative 
methods, such as the direct measurement taken by the SA 
inflight. 

To investigate the correlation between the motion patterns of 
light spots and the collimator state, we analyze three scenarios. 

Case 1: The movement of the three frosted glass light spots 
is influenced by two factors: first, the relative rotation a and 
translation D that occur between the two grating substrates, and 
second, the axis offset 8 of the DM optical axis, which 
corresponds to the schematics shown in Figure 3(a). 
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Figure 1. Layout of SAS in HXI. 
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Figure 2. The working principle of SAS. 


Case 2: The motion of the DM solar light spot is due to two 
factors: the angular shift y occurring in the SA optical axis 
(collimator main axis); or the offset 8 occurring in the DM 
optical axis, as shown in Figure 3(b). 

Case 3: The movement of the SA solar light spot is 
influenced by two factors: the angular shift y in the SA optical 
axis (collimator main axis) and the misalignment T in the SA 
optical lens concerning the detector, as shown in Figure 3(c). 

To enhance our ability to create precise inversion equations, 
it is imperative that we calibrate the optical parameters of the 
DM and SA subsystems. Unlike the parameter calibration of a 
traditional vision measurement system, our focus is on the 
relationship between the motion mode of the target and its 


corresponding light spots to satisfy the SAS function. It is 
proper to directly calibrate the intermediate parameters 
associated with the upper four variables in addition to the 
basic internal and external parameters, principal spot position, 
and aberration. The primary intermediate parameters are: 


1. The proportionality K of the image point’s movement 
resulting from the target’s movement on the front 
substrate (such as hair glass or SA lens), as illustrated in 
Figure 4. Considering that the SAS system contains three 
frosted glasses, allowing for the utilization of the average 
of the three sets of quantities K = (7h;/H) / 3. 

2. The relationship between the angle shift of the DM optical 
axis and the corresponding movement of the light spot. For 
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Figure 3. Three types of light spots moving mode of SAS. 
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Figure 4. Parameter calibration with three frosted glasses. 


the Sun target, it denotes by Lp, and has Lp = xp/tan 0p; 
For the frosted glasses target, it denotes by Lg and has 
Lg = dl/d0g = Lp - sec? 9g, as shown in Figure 5. 

. The relationship Ls between the angle shift of the SA 
optical axis and the corresponding movement of the light 
spot, which can be calculated by Ls = xs/tan Os, as 
Figure 6 described. 


The four intermediate parameters will constitute a crucial 
component in the development of SAS inversion equations. 
Further elaboration on these details will be provided in 
Section 3.3. 


3. Data Processing and Regularization 


Combining the physical model of the SAS, an inverse 
algorithm has been designed to compute the upper four 
variables. Due to the calculation complexity being beyond the 
capacity of the FPGA, the algorithm does not operate in flight. 
We just need to record the light spots’ positions in orbit as 
input data and correspond it to some certain time point, such as 
the flare eruption, and then the inverse equation is solved on the 
ground subsequently. Thus, we obtain the instantaneous states 
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of the HXI collimator each time we are concerned. In the 
following context, the specific inverse equation set will be 
demonstrated and the ideology of how to resolve the equation 
will also be introduced. 

In practice, the algorithm requires two types of input data: 
the position benchmark of the light spots’ original location and 
their new position which occurs at a subsequent moment when 
changes may happen in the HXI system. Although five light 
spots can construct five formulas, only four unknown quantities 
need to be calculated, resulting in an overdetermined problem. 
Despite the possibility of solving it with a common numerical 
method, its ill-posed feature deems them all unfeasible. For 
this, the inverse process is divided into two subprocedures. The 
first involves computing the rotation and a so-called coupled 
displacement of the two substrates using the three frosted glass 
positions. The second involves determining the remaining 
variables using the outcomes from the first subprocedure and 
the positions of the solar light spot. Further explanation of these 
two subprocesses follows. 


3.1. The Over-determined Function 


The first subprocess that calculates the state between the two 
grid substrates is realized by the DM. Obviously, it is a 
common SLAM model that projects the object coordinate to 
the camera coordinate, the equation can be demonstrated as 
follows: 

pi=K-p. (1) 
In the approximate case of neglecting the Z direction, p = [x, y, 
1]” represents the reference position of the three frosted glasses, 


while p’ = [x’, y’, 1]’ represents the deformation occurrence, 
they are both the input data. The K demonstrate the DM 
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Figure 5. DM parameter calibration. 


projection matrix with the form of: 


cosa@ —Sin@ Cy 
K = | sina cosa cy (2) 
0 0 1 


where the quantity a represents the rotation of the two grid 
substrates, and C = [c,, cA is a so-called coupled displace- 
ment coming from D and 8. Combine the three light spots, 
Equation (1) comprises six formulas with the unknown 
quantities a, Cx, and c,, it constitutes an over-determined 
question. Compared to the full rank equation, the over- 
determined problem here cannot achieve the absolute accuracy 
solutions but still has the advantage of obtaining a moderate 
solution that adapts to the measurement errors, such as the least 
squares (LS). 

However, if the LS method is used to solve for (1) to 
evaluate the relative motion between grating substrates, it is 
necessary to determine the center point of rotation and 
translation first, otherwise different definitions of the center 
position will yield different translations, even to the extent that 
at a particular center of rotation the translation is 0, which 
corresponds to the front grating substrate sweeping around that 
point by an angle. 

The definition of the center of rotation on the front grating 
substrate is based on the theoretical center of the DM attenuator 
window. However, as this position cannot be measured 
directly, a substitute is used, namely the center of the line 
between the two frosted glasses, denoted as p, in Figure 7. 
Consequently, the rotation of the front grating substrate 
corresponds to the angle of rotation of all points around pe, 
and the translation corresponds to the amount of movement of 
the pe point. Specifically, the objective function in (3) is 
directly used to solve the rotation based on the location of the 


Ls 


Figure 6. SA parameter calibration. 


three frosted glasses, we have: 


Q = argmin||p’ — Kpl, i=1,2,3 (3) 


where p; is the initial position of the three frosted glasses, and 
pi is the new position of the three frosted glasses. Then the 
coupled displacement is defined as: 


E <p W 
where pe = (pı + p2)/2, and p! = (pi + pa)/2 


3.2. The Rigid Transformation Model 


In actual use, we discover the light spot position measurement 
errors that come from the center extract algorithm, or the 
instrument systemic errors, are a random quantity, which cannot 
be predicted and corrected. It introduces a non-existent deforma- 
tion among the light spots’ graph, which would mislead the LS 
solving process and cause an unexpected wrong value sometimes. 
The problem can be described as a non-rigid transform. To tackle 
the problem, we cite an ARAP iterative method that calculates the 
rigidity of the local transformations from the surface domain. The 
calculation results of this method can also be used as a 
comparative validation of the LS method. The ARAP separates 
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Figure 7. Patterns of rotation and translation of the front substrate. The rotation centern (green dot) is defined as the center of the connecting line of the frosted glasses 


1 and 2. 


Figure 8. The cell structure on the front substrate. The red spots are the three frosted glasses, and the green spot is a virtual point that represent the rotation center. 


the object graph into small cells, and combines the vertexes of 
each cell forming the “energy of cell,” and the energy is 
demonstrated as Equation (5), 


E(Ri, p) = X will — p) — Ri; 


JEN (i) 


PPG = 1 M) 


(5) 


where M is the sum of the points in the cell, p; = [x;, y,|’ is the 
known initial position of ith point, p; =e; yT is the new 
position after transformation, R; denotes the rotation matrix that 
happened through the ith point in the model, j € N(i) represents the 
jth point that connecting to ith point, and || - || denotes the general 
norm. We can use Equation (5) to compute each cell point rotation 
R =[R;] and the coordinates p’ = [ P| after the transform, the 
more minor the energy E, the higher rigidity of the cell shape. 
Constraining the cell energy minimum ensures that the transforma- 
tion is computed with the cell shape as rigid as possible. It would 
eliminate the non-rigid effect to some extent and the solution 
accuracy and robustness would be increased considerably. 

In the SAS, the frosted glasses form a triangle cell on the front 
grating substrate surface, as shown in Figure 8. The frosted glasses 
spots’ positions p; and p; are known input data, which means they 
could be regarded as constrain points. What is more, a virtual 


“controlled point’ that moves following the “constrain points” 
passively is desired to be defined as the rotation center, but also a 
translation center of the cell. Here we choose the geometry center of 
the frosted glasses pı and p2 as the controlled point p, as it simply 
corresponds to the real center of the front substrate. Finally, we 
obtain a cell with four points whose variations can explicitly 
characterize the deviation between the two grating substrates. The 
deviation can be divided into two parts, first is the rotation through 
the controlled point, and the other one is the displacement that 
corresponds to the movement of the controlled point. 

Equation (5) is to compute two variables R and P; to 
minimize E, it is a nonlinear optimization problem that hard to 
find the global optimum. For simplicity, we utilize a Quasi- 
Newton method that fixes one variable value as known and 
calculates the other, then makes an exchange to realize the 
iteration to obtain the global optimization gradually. In 
practice, we first set an initial value for the controlled point 
and combine the observation data to form the p’, then the cell 
vertex vectors can be defined as Equations (6) and (7): 


“= Pi —B, ©) 


ej = p! — pi. (7) 
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By regarding the points’ position values p’ as known 
quantities, Equation (5) can change into matrix form: 


M 


1 T x 
E(R;) z= > wi (ej = Riej) (e; = Rie;j) 
i=1jEN(i) 
u T T T 
IT 1 1 T 
=N wyle; ej — 2e; Riej + ej ej) . (8) 
i=1jEN(i) 


w12 + w13 + w14 —Ww]2 

Ac —w2] W21 + w23 + w24 
—w3] — W32 
— w41 — w42 


As Equation (8) is the function of variable R;, the minimization 
of E only concerns the term R;, which is 


E(R) = 5 (-2wijej" Rie) 
JEN (i) 


= remak Z (oe) 


JEN (i) 
= argmaxTr(R;S;). (9) 


The e; and ej are column vectors that build a weighted square 
matrix S, where the 


IT T 
Si = x Wij Cij Cj = U; D; V; 
JEN 


and the matrix U;, D;, and V; are derived from the singular 
value decomposition (SVD). In Equation (9), according to the 
rule of matrix trace, R; would be chosen to satisfy the 
maximum of Tr(R;S;) when keeps the R;S; a symmetric positive 
semi-definite matrix, thus, here the R; value has: 


R; = VU. (10) 


It is the first iterative step we applied to obtain the rotation 
matrix R, and the second step is to use the known R value to 
deduce the new cell vertexes positions p’. At this time R is 
fixed, we get the p’ by maximum the E value by the partial 
derivative calculation: 


JE (R, p’) 


a = 0. (11) 


Equation (11) is a linear equation set contains M formula, 
which the ith formula has the form: 


D ve -p= Vo 


JEN (i) JEN (i) 


ZUR; + R); — p). 02 
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While the unknown quantities p! in Equation (12) form a 
column vector, it can be transformed into a matrix form: 


Ap' =b (13) 


where the coefficient matrix A elements are weighted by w;j: 


—Ww13 W114 

—W3 —W24 (14) 
w31 + W32 + W34 —W34 

— w43 War + W42 + w43 


and b is the right side of Equation (12). Actually, we only need 
to compute the controlled point position pl, as the constrained 
points position Py (k € constrain) is already determined by the 
input data. Thus the corresponding kth row and column of A, 
combining the kth row of b, are required to be removed, only 
leaving the elements concerning to the controlled point, then a 
new equation is reconstructed as: 


Ap! = b. (15) 


By solving Equation (15), the unknown controlled points p! 
are computed to satisfy the lowest energy of cells. In practice, 
we have three constrained points corresponding to the frosted 
glasses, as a result, Equation (15) has degenerated into a single 
element equation and the unique variable’s solution p! 
demonstrates the new position of the rotation center. 

Finally, the new rotation matrix R and new cell vertexes 
p= [p, P p] are brought back into Equations (15) and (10) 
separately to perform the iteration process and stop at a suitable 
condition, such as when the amount of iterative change is less 
than a certain threshold. We use the rotation center to describe 
the rotation and translation of the substrate, thus, R = R, = a 
and C = p! — p. 


3.3. The Inverse Function 


As previously stated, the ARAP algorithm provides solely a 
relative rotation a and coupled translation C between the front 
and rear substrates. To solve for the remaining state variables 
D, G, and y, it is essential to establish the equations along with 
the motion pattern. Examining the first motion pattern in 
Figure 3(a), we can see that both D and ĝ cause coupled 
translation C, enabling us to obtain the first equation of motion, 


D+ß-Le=C=p -—p,. (16) 


Similarly, according to the second mode of motion 
Figure 3(b), the motion of the Sun light spot in the DM is 
jointly caused by 8 and y, and the second equation of motion 
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can be obtained: 
(6+ 7)-lp=q'-@ (17) 


where q is the initial position of the DM Sun light spot and q’ is 
the position after the instrument state change. According to the 
third motion model Figure 3(c), the motion of the SA Sun light 
spot is also synthesized by two factors, y and T, respectively. 


T+y L=s'—s (18) 


where s is the initial position of the SA solar light spot, and s’ is 
the position of the spot after a state change occurred in the 
collimator. T is a special variable representing the offset 
produced by the lens in SA placed on the front substrate 
relative to the detector mounted on the rear substrate. Since the 
SA lens is also a target spot in the DM and is itself a quantity 
measured by the DM, T can be expressed in terms of the 
rotation and translation of the front substrate: 


ti —+¢ 
B 


where ¢ is the initial coordinate of the SA lens captured through 
the DM, and f’ is its new coordinate after translation and 
rotation followed by the front substrate, which has 
t = Rt + D. The negative sign indicates the object-image 
opposite relationship. Note that the tangent calculations are 
simplified to arc angle calculation in Equations (16)-(18), 
taking into account that the actual deformation of the collimator 
is very small. Combine Equations (16)—-(18), we determined the 
measurement equation set of the SAS: 


By =g (20) 


T=- 


(19) 


where 


1 Lg 0 
B= 0 Lp Lp 
—1/K 0 Ls 
and its elements involve the optical parameters defined in 


Section 2. y = [D (34]’ is the vector to be solved. The right side 
vector 


oe | 
R-t-t 
K 


gs = 
si — s+ 


represents the observed data. So far, by back-calculating 
Equation (20), 
y=B'g (21) 
we realize the calculation of the state variables D, 8, and +. 
3.3.1. Tikhonov Regularization Method 


Equation (21) is a full-rank linear equation, commonly it can 
be straightly solved by Gaussian elimination or the Gaussian 


Yu et al. 


Seidel method. Unfortunately, inspecting the condition number 
of matrix B described as cond(B) = ||B™'||, - ||B|l2, which 
illustrate the illness of a linear equation. A larger condition 
number represents a more unstable inverse solution to the 
equation and a more sensitive response to noise, such as, in 
Equation (20), we get an enormous condition number bigger 
than 440,000. To analyze the reasons, we research the solution 
structure with the SVD method, by matrix operation, the 
solution of the equation can be given by 


T 
y=B g=) En, (22) 
i Gi 

where u; and v; are the orthogonal column vectors of the SVD 
left and right matrix of B. g; is the singular value. It is obvious 
that the coefficients u ¢/o; are impacted extremely by ø; as it 
would decay into a tiny value in the ill-posed matrix. We 
impossibly get an ideal input g with no noise but a data g + Ag 
with a perturbation, as a consequence, in the coefficient 
uf (g + Ag)/o;, the perturbations Ag is magnified and causes 
the theoretical exact solution y becomes “unbounded.” Thus, 
the traditional method cannot be used in the ill-posed problem 
computation, but some specialized approach is required. 
Equation (20) can be expressed as a least squares (LS) form: 


Q = min{By — g}. (23) 


To restrict the range of possible solution values, a penalization 
term ||y|| is introduced in Equation (23), which is known as 
regularization. Then Equation (23) becomes the following 
form: 


Q = argmin {||By ~ gl; + Xll} (24) 


Consequently, the solution of Equation (24) can be reformu- 
lated as 


T 
U: 
yy = (B7B + XIB" = D g'n (25) 
i i 


where gl^ = o? / (o? + X). It adds a controlling weight that 
satisfies the demand, sometimes the method is called L-2 
regularization or Tikhonov regularization. The A is a positive 
parameter that balances the two parts, the residual norm 
\|By — g||2 and the solution norm ||y||2. The weight is enhanced 
to diminish the solution norm |ly||2 by enlarging the A value, 
which can limit the uncontrollable perturbations in 
Equation (22) to some extent, but may cause a bigger deviation 
of between the regularization solution and the ideal solution 
due to the limitation of ||By—gll2 is expanded. On the 
contrary, if the A is reduced, the regularization effect is 
weakened. The choice of À value then becomes the core of the 
problem. There are many well-established methods for the 
selection of regularization parameters, such as the Discrepancy 
Principle (Engl 1987), the L-curve Criterion (Hansen 1992; 
Johnston & Gulrajani 2000), the Generalized Cross Validation 
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(Golub & von Matt 1997), and the Normalized Cumulative 
Periodogram (Mojabi & LoVetri 2008), etc. In this paper to 
compute the appropriate À, we utilize the discrepancy principle, 
according to Equation (25), we can obtain: 


2 T 
Oj uj g 
By, — g =USV' . : Levi g 
i 2 te G i 
X2 
= SE ai e (26) 


Then the residual norm is: 


x i 
lB», — gle = o i us| . (27) 
Taking into account the observation noise, there is 
g = g™°" + e. Based on the properties of SVD decomposition, 
the coefficient u7 g°***t follows a monotonically decreasing 
pattern and satisfies the discrete Picard condition. There exists a 
bounding value t, when o;< t, because of the dominance of 
uj gt, the uf g = uf (®t + e) continues to decrease 
monotonically. However, the noise vector e may not align with 
gt, which means that uï e does not necessarily satisfy the 
discrete Picard condition. Consequently, when g; > t, the noise 
e is larger than g®*“', then the vector u e will dominate leading 
to a platform or bounce in the value of u;’g instead of 
decreasing, which would magnify the superimposed solution y, 
with tiny singular values. Hence, the main aim of the 
regularization process is to reduce this adverse effect, and the 
parameter y! Al works as the modifying factor. To analyze the 
effect of regularization parameters, we assume that the 
observation noise e is the Gaussian distributed with 


Cov(UTe) = UTE(eTe)U = fI, 


7 is the standard deviation of e, then Equation (27) can be 
expressed as 


By, — g = [a - put gf 


dn 
2 
SPA- put ge + @ AA — ng 
i=1 


(28) 


where A, is the boundary of the platform region and we have 
E(\lel|3) = ||n||5 that represents the expected value of the noise 
vector. As Equation (28) expression, the residual norm is 
composed of two parts, representing the effect of the 
observation error and the exact solution value on it, 
respectively. The ||By, — g| increases with the À value 
increases, and e has a smaller value than g* commonly, 
whose proportion in the residual norm gradually diminishes 
along with the increasing of À. Then we analyze two scenarios: 

Case 1: there exists a critical value 7, when A<7, the 
singular values slightly increase, and the divergence of the 
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superimposed solution u; g/o; remains dominated by noise e. 
This indicates under-regularization (under smoothing) and the 
solution y = yr", 

Case 2: when \ > 7, g; increases significantly, which largely 
suppresses the amplitude of the superimposed solution u’ g/o;. 
It demonstrates the over-regularization and the residual norm is 
collectively determined by both observation noise and exact 
observation value. 

The discrepancy principle is to find a critical A such that the 
residual norm is exactly dominated by the observation error e, 
ie., there is ||By, — gl} = |le|3, which determines the 
regularization parameter that achieves equilibrium between 
the residual norm and the solution norm. This article advocates 
the use of the discrepancy principle when selecting the 
regularization parameter, owing to the ability to pre-evaluate 
the observation noise e (equivalent to the spot center extraction 
accuracy). Nevertheless, this still exists a potential risk, as an 
inadequate estimation of the observation error could signifi- 
cantly impact the solution outcomes, particularly when the 
parameter estimation is low, resulting in under-regularization 
that amplifies the range of the solution norm and causes the 
deviation of y^ from the ideal solution y to be excessively 
large. To ensure caution, a safety factor v around 2 is typically 
added to the noise for a slightly over-regularized state, then we 
have: 


By, — gle = v linll- (29) 


According to Equation (29), we compute and substitute it 
into Equation (25) to obtain the regularized solution. 


3.3.2. In Bayesian Perspective 


Expanding on the Tikhonov regularization method discussed 
in Section 3.3.1, the Bayesian approach (Jin & Zou 2008a) 
offers fresh insight and makes a cross-validation of the solution 
of Equation (23). Assuming a prior distribution for the 
observation noise and the target variable, one can estimate 
unknown information by determining the posterior probability 
density function (PPDF). In the probabilistic perspective, the 
solution aims to determine the probability value that maximizes 
y, taking into account the known observation g, and uses the 
maximum a posteriori (MAP) of the y to estimate the optimal 
value. Based on the basic arithmetic of Bayesian estimation, 
there is 


m(gly)7(y) (30) 


m(ylg) = ae) 


where m(g|y) is the likelihood function, z(y) is the prior 
probability of y (Jaynes 1968), and 


m™(g) = frero 


is referred to as a normalized constant, which is usually ignored 
in calculations, then Equation (30) can be expressed as a 
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proportional relationship: 


Tyle) x m(gly)7(y) (31) 
thus, MAP can estimate the y value: 
Ymap = argmax {7(ylg)}- (32) 


To resolve the issue at hand, it is essential to acquire the 
probability prior and the probability likelihood as denoted in 
Equation (31). Initially, considering the noise of the observa- 


tion, Equation (20) can be expressed as 
g= By +e. (33) 


The observation noise e follows a continuous probability 
distribution at R”. Assuming the independence of e and y, the 
observation g is consistent with the noise distribution 


(34) 


Usually, the noise e is the independent identically distributed 
Gaussian random variables with mean zero and standard 
deviation gı. Equation (34) is given by 


20? 
|By — sÈ) 


20? 


m(gly) = m(e) = T (g — By). 


1 
T(ely)= vas 
1 


Ta a| (35) 
Defining the a priori probability m(y) presents difficulty as the 
quantity of deformation y in SAS measurements is unknown. 
Consequently, we should assume the HXI instrument tends to 
stabilize, with minimal deformation in the absence of other 
interferences, and a lower likelihood of larger deformations. It 
thus follows a Gaussian distribution with parameter (0, o2) 


1 y? 
T) = exp| -= |. (36) 
V2T 02 o( 5) 
Combining Equations (31), (35), and (36), we obtain 
By — glP llylP 
T x ex ex 37 
Olg) il 20? p 207 (87) 
then the object function of MAP is 
By — 2 2 
Q = arg max T (y|g) = arg max exp |By sl Ixi 
y y 20i 205 
(38) 
or another form 
_ f ilBy — gl? | lbi 
= arg min {| ———~——. + —— >. 39 
= j y í Qa? 203 


Comparing Equations (24) and (39), it reveals that the 
regularization parameter resembles the probabilistic model 
parameter utilized in the Bayesian framework. We introduce a 
hyper-parameter in Equation (31), which is also a random 
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variable and helps to establish the correlation between the 
Bayesian model and the classical Tikhonov regularization 
theory. According to the Bayesian theorem, the PPDF is 


m(y, Alg) x m(gly, A)» TVA) > TO). 


As the observation g is not affected by the hyper-parameter A, 
the first term on the right-hand side has z(gly, A) = (gly). If 
the coefficient matrix B € R”*", considering the interrelation- 
ships between observations, there exist 


(40) 


g)" S-\(By — 9} 


(41) 


m(gly) = 


er 
(27)? [det (S) 2 


In this article, g is composed of the computational center of 
three types of light spots, and the observation noise reflects the 
extraction accuracy of light spots. Since the components are 
independent, the noise has a covariance matrices as 


Ne 
S= Na A 
ns 


where ne, Na and 7, are the error variances of the three vectors 
of the observation data. Due to the regularization parameter 
being a penalty factor applied to the solution norm ||y||2, it is 
manifested as a correction of the prior probability z(y|A) within 
the Bayesian framework. As y follows a Gaussian distribution 
and its components satisfy the covariance matrix W (Calvetti & 
Somersalo 2018), for whom the A serves as a scaling factor and 
ultimately modifying the solution norm of the y. In this 
scenario, y satisfies the Gaussian prior 7(y|\) ~ N(O, AW) with 
a hyper-parameters 

(yA) = S exp( -7w (42) 
(27)? - ydet (AW) 2A f 
considering that the components of y are independent of each 
other, then W must be a diagonal matrix. It is important to note 
that this paper analyses the different deformations represented 
by D, 8, and y according to the design specifications of the 
SAS, combined with certain empirical assumptions and 
simulation data. The argument is made that they should 
conform to different Gaussian priors. This involves estimating 
the variance of the deformation distributions based on the setup 
of the simulation inputs. It is stated that 


where oy, 73, and o, represent the standard deviation of the 
three preset deformations, respectively. m(A) represents the 
prior distribution of the hyper-parameters and can be described 


by the conjugate prior (Gelman et al. 1995). For this study, we 
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Figure 9. Different types of light spots in SAS in the ground test. Some simulated light sources are desired in the ground test. For frosted glasses, we used three LED 
light sources mounted only before the frosted glasses; a light pipe is the alternative to the incident solar light, and the quite tiny light spots are caused by its small FOV. 


will use the gamma distribution with parameters a and 8 


T) = LxA 


4 
Tra) (43) 


where I is the gamma function. Bringing Equations (41), (42), 
and (43) into Equation (40), we obtain the joint PPDF of y and 
À 


1 
m(y, Alg) x exp(—5(B — g)'S\(By — 9) 


T an 


1 
x4 apf- 2) “A Texp(—GA) (44) 


then the objective function of MAP estimation can be easily 
determined 


Tyy7-1 
Q = argmin | +| cay — gyts-\By — gy) + YY 
yA (2 À 


ce E i= a) ina $ a} (45) 
The Equation (45) comprises two optimization variables, y and 
A, which can be resolved via an iterative approach (Jiang et al. 
2020). 

Step 1: Define the initial value of the parameter Ao, and 
subsequently, y for the kth iteration is evaluated in accordance 
with according to Equation (25) 

w-! 


-1 
——]| BS~'g. 
% ) i 


y= Ga 4 (46) 


Step 2: Update the parameter ,,, by differential operation 
0Q/OAly=y, = 0, then Ax+1 satisfies the function 
ye wW- p7 
2Ai+1 


age eet 
+ B= 0. 


(47) 
Akt 


The iteration ends when the abort condition ||y, — y,—1|| < € 
is met. We now achieve the estimated solution y with the 
Bayesian model, which is an approach that differs from the 
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Tikhonov regularization and denotes a cross-check of the 
certainty of the solution. 


4. Numerical Simulation 


For the procedures discussed previously, it is necessary to 
conduct simulations to assess their effectiveness in regularizing 
the ill-posed inverse problem in SAS. For this purpose, 
simulated data will be used to validate the accuracy of the 
inversion, and the results of real tests will be presented later. 
However, before this, we first calculate the light spot extraction 
accuracy of SAS, which serves as the observation error. 


4.1. Accuracy of Light Spot Position 


During both ground testing and inflight operations, SAS 
detects various types of light spots due to different illumination 
sources, as demonstrated in Figures 9 and 10. Different center 
coordinate extraction algorithms are implemented for the 
various forms of light spots, based on their distinguishing 
features. 


1. Circular light spots with symmetric shape and uniform 
energy: binarized center of mass algorithm; 

2. Bessel-shaped circular light spot: ellipse fitting center 
method (Fitzgibbon et al. 1999; Zhou et al. 2018); 

3. Full-disk large circular light spot: four-quadrant area 
fitting algorithm or 14-row edge data fitting algorithm; 


The principles behind these algorithms and the accuracy of 
the extraction will be thoroughly described in our forthcoming 
publication. Due to space constraints, we will be citing the 
findings and presenting the results in Table 1. It is evident that 
the extracting precision is equal to that of the noise level in the 
observed data g, which also serves as a crucial input for 
subsequent simulation calculations. 


4.2. Operation Result 


The numerical simulations conducted in this paper were 
executed on a computer equipped with an Intel Core i7 CPU 
and 16GB of RAM, utilizing Python 3.8.5 to execute the 
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Figure 10. Different types of light spots in SAS inflight. 


Table 1 
Characteristics of Types of the Light Spots in SAS 


Light Spot in SAS Scenario Morphological Structures Extracting Method Accuracy (Std) 

Frosted glasses Inflight use Round spot Shape centroid 0.04 pixels 
Ground test Round spot Shape centroid 0.04 pixels 

Solar image in DM Inflight use Round spot Shape centroid 0.04 pixels 


Ground test Multi annulus 


Circle disk 
Multi annulus 


Inflight use 
Ground test 


Solar image in SA 


Ellipse fitting 0.025 pixels 


Four quadrant method or 14-row fitting 
Shape centroid 


0.11 pixels 
0.04 pixels 


algorithms. For all subsequent simulations, we established a 
consistent initial position go, for the spot and generated 1000 
sets of deformations randomly to be implemented as simula- 
tions. The sequence employed for the simulations is as follows: 

Step 1: Utilizing the imaging principle of SAS, we first 
derive the altered position gı of the light spot caused by 
simulated deformation. The resulting value g=g)— go, 
represents an ideal observation free from noise. It should be 
announced that the rotation and translation center of the front 
grating substrate is defined as the center of the line joining 
frosted glass 1 and 2 during the forward derivation of the 
changing coordinate g4. 

Step 2: To the ideal observation, we add an observed error 
gZ = g + e, where e represents the extracting accuracy of the 
light spot center. 

Step 3: Solving the rotation angle and coupled displacement 
with LS or ARAP method. 

Step 4: Solving the inverse problem Equation (21) by 
utilizing the Tikhonov and Bayesian methods and subsequently 
evaluate the accuracy of the inversion algorithm. 

It should be noted that we have made an assumption that any 
deformation in the HXI collimator during its full life cycle will 
not surpass the permissible threshold that would affect the 
quality of the flare imaging when conducting the numerical 
simulation. This constrains us to set the peak of simulated 
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deformation equal to the maximum allowable deformation of 
the collimator, whereby a=+ 10” (the threshold that could 
impact imaging details), d= + 18 um (matching the period of 
the thinnest grating slit), 8=+10”, and y=+ 10". 

For comparison, it is assumed that there is no noise (i.e., 
e =0), the results obtained from directly solving Equation (1) 
with the LS and ARAP methods, and Equation (21) using the 
Gaussian Elimination (GE) method, are presented in Figures 1 1 
and 12. At this condition, the result demonstrates the high 
accuracy that reaches almost the zero inversion error. The LS 
method has higher inversion accuracy compared to the ARAP 
method, which is due to the fact that the ARAP does not solve 
the equations exactly, but instead calculates them through 
iterative optimization. These high accuracies show the correct- 
ness of the inverse process modeling. 

Figures 13 and 14 display the outcomes upon adding the 
observation noise e into observation data g. The two methods, 
LS and ARAP, have very similar arithmetic results. The peak 
error of œ is about +15”, and the coupled displacement is about 
0.1 pixels (as the C is an intermediate process variables, we 
demonstrate it in the “pixel” unit). Similarly, we analyze the 
error values at the 1/4, 1/2, and 3/4 quantile, which are 
indicated by each of the six black dashed lines in 
Figures 11-19. At the 1/2 quantile, the error of a is about 
+3"5, and the error of cą and cy is about 0.03 pixels. In the 
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Figure 11. Inversion error distribution with the LS and GE methods without observed noise. 


following numerical simulations, we use only the ARAP 
method to calculate the rotation angle for simplicity. 

As a consequence, due to the ill-posed nature of the inverse 
function (21), an unregularized computation of the D, 3, and y 
variables, such as GE, results in a highly deviated and 
unacceptable solution that is far from ideal, as shown in 
Figure 15. 

When using the Tikhonov method to solve Equation (21), 
the results are displayed in Figure 16. The a calculation is 
unaffected and retains the same accuracy of about 3”5, and 
regularization significantly optimizes the solution outcomes for 
D, ß, and y, whose peak error is approximately +55 um, +10”, 
and +10”, respectively, and the errors displaying a uniform 
distribution within their peak range. At the 1/2 quantile, the 
solution error is estimated to be around D = + 24 um, 

=E aie and {= 25 5! 

The solution results utilizing the Bayesian framework exhibit 
distinct characteristics. The PPDF, which takes Gaussian noise 
and Gaussian prior distribution into account, generates an 
inverse error with a Gaussian distribution structure of a mean 
value of 0. This suggests that computation based on the 
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Bayesian framework is more likely to achieve smaller inversion 

errors, as illustrated in Figure 17. Based on the findings, the 

MAP produces slightly more significant peak errors than the 

Tikhonov approach. Specifically, the peak values for D is about 

+65 um, 3 and y are both slightly exceed +10”. However, the 

error values at the 1/2 quantile are lower, with D = +20 um, 
= +3”, and Vi +3”, 

We have observed that the maximum error of the inversion 
changes with the preset simulated deformations. Specifically, 
when the preset deformations are large, the final inversion error 
increases, whereas it decreases when the deformations are 
small. Besides, the distributional characteristics of the error 
persist without alteration. For instance, the conclusions 
presented in Figures 16 and 17 are the outcomes acquired 
from a predetermined deformation of a within +10”, D within 
+18 um, 6 within +10”, and y within +10”. By applying 
another set of predetermined values of a within +5”, D within 
+5 um, 6 within +3”, and y within +3”, the inversion 
outcomes are exhibited in Figures 18 and 19. Compared to 
larger preset deformations, the computation of rotation does not 
present an ill-posed problem. Therefore, solution accuracy 
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Figure 12. Inversion error distribution with the ARAP and GE methods without observed noise. 
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Figure 13. Rotation and coupled displacement error distribution using LS. The equation is solved with the Powell and BFGS methods, with the same result. 
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Figure 14. Rotation and coupled displacement error distribution using ARAP. 
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Figure 15. Direct inversion error distribution with the GE method adding observed noise. 
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Figure 16. Inversion error distribution with the Tikhonov method. 


remains unchanged. The peak error of D, 3, and y solutions is 
lower with the implementation of the Tikhonov method, which 
achieves a peak error of about D = +20 um, 8 = 3”, and 

=+3”, and a 1/2-quantile error of D=+8pum and 
B= 1”5, y= 1”5; The MAP method can achieve the highest 
of about D = +25 um, B= +4", y= +4", and at the 
quantile, the error values are D= +6 pum, G=+1", 
y= +1”. The error distributions of both methods are unaffected 
by the magnitude of the preset deformation. The reason is that, 
when carrying out regularization, the computation incurs a 
penalty factor in the form of the solution norm ||y*|ĝ, and for 
the identical inverse equation, the ideal solution norm increases 
as the preset deformation increases. At this time, to meet the 
optimization objective of Equation (24), the solution norm 
contour expands and intersects the smaller residual norm 
contour—resulting in a smaller ||By — g|}. Such regularization 
is inclined toward under-smoothing, leading to the inverse 
problem being closer to an unbiased solution. Hence, it 
achieves the “divergent” solution that has larger amplitude. In 
summary, a significant deviation in the collimator would 
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decrease the precision of the SAS monitoring, which 
necessitates the assumption that the instrument would not 
exhibit a substantial deviation, as indicated at the beginning. 
Furthermore, if the collimator indeed undergoes an unaccep- 
tably large deformation, which will compromise the imaging 
quality of the HXI instrument, it may not be necessary to 
prioritize the measurement accuracy of the SAS if the usability 
of the HXI cannot be ensured. Fortunately, the HXI underwent 
precise calibration using an X-ray beam, additionally, the 
installation and integration processes were in a relatively stable 
circumstance, and the instrument was accurately temperature- 
controlled inflight to ensure an unchanged state. Based on these 
factors, it is reasonable to assume that significant deformation 
would not occur in the collimator. These assumptions and 
conclusions were subsequently confirmed through inflight 
testing. 

In summary, the simulation results are presented in a 
tabulated form for clarity. Table 2 compares the results 
obtained from inverting simulations using both Tikhonov 
regularization and Bayesian frameworks for large preset 
deformations. The MAP approach offers some advantages 
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Figure 17. Inversion error distribution with the MAP method. 
Table 2 
Inversion Error with Large Preset Deformations 
Tikhonov Bayesian 
1/4 1/2 3/4 Peak Std 1/4 1/2 3/4 Peak Std 
Quantile Quantile Quantile Quantile Quantile Quantile 
a(arcsec) 1.74 3.57 6.10 16.57 5.22 1.74 3.57 6.10 16.57 5.22 
dum) 11.68 23.86 37.46 55.42 28.89 10.09 19.66 31.14 65.90 25.46 
dum) 11.85 22.92 38.22 57.62 29.11 9.33 19.60 30.64 63.87 25.34 
3, (arcsec) 2.39 5.11 7.60 9.99 5.84 1.68 3.32 5.32 12.29 4.41 
By (arcsec) 2.86 5.06 7.70 9.99 5.90 1.55 3.24 5.21 11.53 4.37 
7, (arcsec) 2.53 4.85 7.35 9.99 5.68 1.70 3.29 5.28 11.09 4.32 
7, (arcsec) 2.43 4.80 7.40 9.98 5.71 1.58 3.27 5.23 10.83 4.31 


owing to the Gaussian distribution which guarantees a higher 
proportion of small errors, the results of these two methods will 
be jointly considered. Table 3 presents the inversion outcomes 
obtained with smaller preset simulated deformations, some 
corresponding values are smaller compared to Table 2. Since 
the simulation data utilizes 1000 groups of randomly preset 
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values, each simulation could produce slightly varied out- 
comes, but the results remain at the same level. It is worth 
noting that the discrepancy principle is capable of reducing the 
peak error when it is over-regularized by increasing the safety 
factor v. However, there is a potential risk of raising v 
arbitrarily, due to the fact that it would constrain the inversion 
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Figure 18. Inversion error distribution with the Tikhonov method upon small preset deformation. 
Table 3 
Inversion Error with Small Preset Deformations 
Tikhonov Bayesian 
1/4 1/2 3/4 Peak Std 1/4 1/2 3/4 Peak Std 
Quantile Quantile Quantile Quantile Quantile Quantile 
a(arcsec) 1.63 3.42 5.82 16.69 5.15 1.63 3.42 5.82 16.69 5.15 
dum) 3.85 7.92 12.32 18.59 9.49 2.96 6.21 10.03 24.78 8.56 
dum) 3.86 7.90 12.27 20.54 957 2.80 6.04 10.61 23.34 8.70 
3, (arcsec) 0.73 1.50 2.14 2.99 1.68 0.52 1.02 1.66 3.86 1.39 
By (arcsec) 0.77 1.48 2.24 3.00 1.72 0.49 1.04 1.70 4.02 1.41 
7, (arcsec) 0.74 1.45 2.20 3.00 1.71 0.48 1.04 1.67 4.07 1.44 
7, (arcsec) 0.70 1.53 2.24 3.00 1.73 0.49 1.01 1.78 4.21 1.47 


solution to a range close to 0. This may lead to the inversion 
error always being the same as the actual deformation, 
rendering the assessment of the inversion error meaningless. 
In addition, numerous practical measurements were con- 
ducted using SAS throughout the HXI assembly, launching, 
and in-orbit operation. Before integration into the satellite, 


collimator, 


third-party measuring equipment was capable of taking on the 
such as the Coordinate Measuring Machine 


(CMM), to validate the SAS comparatively. Figure 20 exhibits 
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the status monitoring outcomes of the HXI at 14 crucial time 
points (illustrate in Table 4). The relatively flat curves in the 
measurement results suggest that the HXI collimator is in a 
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Figure 19. Inversion error distribution with the MAP method upon small preset deformation. 
Tabled less than 6 um in D,, and less than 5 um in D,, thereby 
Measunng Time Pomts demonstrating the effective usability of the SAS. 
State Time Point 
1 HXI assembly completed, initial state of SAS 4.3. Discussion 
2 After HXI vibration test : + ati : : 
The Tikhonov regularization and MAP methods differ in 
3 After HXI thermal cycle test ‘able: hi hei luti l imil d th 
4 After HXI transport test principle, however t er solution results are Simi ar and there 
5 Completion of all HXI tests, from then on the CMM was unavailable must be a close connection between them (Antoni et al. 2023). 
6 Collimator assembled onto satellite Equation (24) shows that the penalty term present in the 
7 Spectrometer assembled onto satellite solution equation is of the form ||y||3, which is usually effective 
: Whole. satellite: assembly, completed but not the only choice. By considering the general form, where 
9 After satellite environmental test th lty t tak hie: f Lyi? H 1989). th 
10 Before satellite launching S penalty erm. a ses the torm l y| (Hansen ), the 
11 Inflight measurement Tikhonov regularization can be expressed as 
12 Inflight measurement 
13 Inflight measurement 
z 2 2 2 
14 Inflight measurement arg E { ||Bx = gll + x I|Zy1| } (48) 


stable state. The SAS solution deviates from the CMM 
measurements, with a less than 5”, D, less than 10 um, and 
D, less than 13 um. Moreover, the two regularization methods 
in SAS present some deviations in results, with an incoherence 
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where L is a discrete approximation of some operation 
operator on the variable y in the first kind of Fredholm 
integral equations (Hansen 1989, 2010), describing the 
interrelationships among the variables to be solved. The 
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Figure 20. Practical use of SAS and CMM in significant time points. The error bars for both the Tikhonov and MAP methods are calculated as 1/2-quantile of the 


10”, D = 


error distribution with the preset deformation of a = 


18 um, 8 = 4 


10”, and y= +10”. Within these 14 measurement points, the initial five points 


occurred before the HXI integration into the satellite platform. This enabled the use of CMM as a comparative measurement. After state 5, HXI can no longer make 
use of third-party measuring devices for assistance, relying solely on the SAS measurement results. 


Tikhonov regularization is a particular case where L=I. 
Upon examination of Equation (45), it is evident that in the 
first term, the residual norm is modified with the observation 
noise matrix S, while the second term outlines the 
covariance matrix W to depict the correlation between the 
observed variables y. Subsequently, if W`! = L7L, the form 
coheres with Equation (48). The form of L is related to the 
definition of the a priori probability, e.g., the Markov 
Random Field used in heat conduction problems (Jin & 
Zou 2008b) defines a priori probability of 
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or similar to the use of a Gaussian distribution with 
hyperparameters in this study, with 


while the y components are independent. In the actual 
measurement, D varies mainly due to the flatness of the 
satellite mounting surface and the collimator deformation 
caused by gravity unloading after launch. Therefore, it is set 
according to the standard deviation of og=18 ym. The 
variation of 3 is primarily due to possible temperature changes 
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in the satellite mounting plane, which cause the front and rear 
grating substrates of the collimator to tilt. This is set according 
to the standard deviation of og= 10". The variation of y is 
mainly due to changes in the observation optical axis caused by 
uncertainty in satellite attitude control. This is set according to 
the standard deviation of ø, = 10”. This implies a constraint on 
the range of y components in different proportions, and 
although this fits the prediction of the a priori probability, it 
does not necessarily achieve the effect of optimal regularization 
in terms of the form of the (48) as the larger the a priori 
variance makes the L smaller, which may result in an under- 
regularization. The assertion is indeed evident from the 
numerical simulation and the actual measuring results. Given 
the above considerations, this research employs the two 
regularization methods that share the common trait of being 
based on the regularization of generalized SVD but differ in 
their defining processes: the Tikhonov method builds upon the 
LS optimization function by incorporating solution range 
constraints, while the MAP method constructs the PPDF to 
build the optimization function through the distribution of 
observation noise and the solution’s priori distribution. Never- 
theless, their optimization functions both strive to balance the 
residuals and solution ranges, thus evidencing a close 
association between the two methods. 


5. Conclusion 


This paper aims at the problem that SAS is susceptible to 
noise from observation data during the monitoring of the HXI 
collimator, resulting in significant deviation of the inversion 
solution. Some methods are proposed to decrease the inversion 
error of SAS by utilizing the ARAP and regularization 
techniques following the physical model of SAS. The 
simulation data demonstrate the validity of the SAS inversion 
model and the efficacy of the regularization process. Addition- 
ally, the results obtained from the real test data provide 
evidence that these methods can realize the online monitoring 
function. In brief, we make some summarizations: 


1. Based on the data simulation, the LS and the ARAP 
methods yield a measurement accuracy of about +15” in 
peak and +3”5 in 1/2-quantile for the rotational 
deformation. The Tikhonov regularization approach 
results in a peak measurement error of about 

= +55 um for the translation, and the peak measure- 
ment errors for the two optical axes are = +10”, and 
y= +10", respectively. Its accuracy at the 1 /2-quantile is 

= +24 um, G=+5", and y=+5”. The measuring 
peak errors from the MAP method are slightly larger than 
those of the Tikhonov method, but the accuracy at the 
1/2-quantile is superior, that is, D = +20 wm, 3= +3", 
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and y= +3”. All the simulated data have considered the 
extraction error of the light spot center; 

2. The accuracy of the SAS algorithm’s inversion is linked 
to its deformation scale. If the deformation is insignif- 
icant, the accuracy of the inversion is high, and 
conversely, if the deformation itself is significant, the 
inversion accuracy is reduced; 

3. Analyzing the several times of practical test data, it has been 
discovered that there are disparities between the outcomes 
of SAS and CMM tests, specifically, a is less than 5”, Dy is 
less than 10 wm, and D, is less than 13 um. Furthermore, 
disparities also exist between the Tikhonov and MAP 
methods, with D, less than 6 um and D, less than 5 jum. 


In summary, this paper demonstrates that the regularization 
method significantly improves the stability of SAS’s inversion 
calculation and enhances its robustness to observation noise 
pollution. Furthermore, the monitoring feature of SAS provides 
vital assistance in completing the HXI mission. It aids in 
instrument assembling and tuning during production and 
monitors the state throughout the entire lifecycle of HXI. 
Additionally, it enables the pointing of the solar orientation 
during in-orbit operations. 

Unfortunately, there are inherent limitations in the visual 
measurement system, for instance, improving the light spot 
extraction accuracy indefinitely is not possible. Also, the 
solution accuracy for inversion equations is determined by 
optical parameters predominantly, making it difficult and 
inefficient to increase accuracy solely relying on data 
processing methods. Therefore, future research should focus 
on designing a more reasonable optical system to improve the 
robustness and anti-interference ability of inversion equations 
for the SAS-like monitoring systems, assisting the data stability 
algorithms to enhance accuracy. The study presented in this 
paper significantly improves the usability of SAS, allowing it to 
successfully perform its intended functions and make further 
progress toward achieving the HXI mission goals. Ultimately, 
we anticipate that such investigations will yield theories and 
techniques that can assist the design of the forthcoming flare 
observation instruments in the future. 
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